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Abstract 



A measurement of the b-hadron production cross section in proton-proton collisions 
at a/s = 7 TeV is presented. The dataset, corresponding to 85 nb^^, was recorded with 
the CMS experiment at the LHC using a low-threshold single-muon trigger. Events 
are selected by the presence of a muon with transverse momentum pj>6 GeV with 
respect to the beam direction and pseudorapidity < 2.1. The transverse mo- 
mentum of the muon with respect to the closest jet discriminates events containing 
b hadrons from background. The inclusive b-hadron production cross section is pre- 
sented as a function of muon transverse momentum and pseudorapidity. The mea- 
sured total cross section in the kinematic acceptance is cr(pp -^h + X^]i + X') = 
1.32 ± 0.01 (stat) ± 0.30(syst) ± 0.15(lumi) ^h. 
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1 Introduction 

Measurements of b-hadron production in proton-proton (pp) collisions at the Large Hadron 
Collider (LHC) are important tests of Quantum Chromodynamics (QCD) in a new kinematical 
region. Results on b-hadron production in proton-antiproton collisions at the lower center-of- 
mass energies, ^/s, of the CERN SppS Collider [IJ and the Tevatron [2-5J have aroused sub- 
stantial interest because of tensions between the experimental results and the theoretical ex- 
pectations in 13- First results at the LHC from pp collisions at ^/s = 7TeV have been reported 
by the LHCb collaboration for the forward rapidity region using semi-inclusive decays [8J and 
by the CMS collaboration in the central rapidity region using fully reconstructed hadron 
decays ||9l. 

The b-quark production cross section in hadron collisions has been computed at next-to-leading 
order (NLO) in perturbative QCD ||T0HT2l . The observed large scale dependence of the NLO 
results is considered to be a symptom of large contributions from higher orders: small-x ef- 
fects pis', ^T^, where x ~ m^, / a/s, are possibly relevant in the low-px domain, while multiple- 
gluon radiation leads to large logarithms of pj/ m-^ and may be important at high px The 
resummed logarithms of px/ next-to-leading-logarithmic accuracy have been matched to 
the fixed-order NLO calculation for massive quarks ftl6fl . At the non-perturbative level, the 
b-hadron px spectrum depends strongly on the parametrization of the fragmentation func- 
tion [17J. The b-quark production cross section has also been studied in the general-mass 
variable-flavor-number scheme ||T8| and the ki factorization QCD approach 119112m . 

In this paper we present an inclusive measurement of the production of b hadrons decaying 
into muons and jets based on 85 nb^^ of data recorded by the CMS experiment using a low- 
threshold single-muon trigger. Muons from b-hadron decays are distinguished from back- 
grounds based on their transverse momentum relative to a nearby jet (p'J^). 

In Section |2] a brief overview of the CMS detector is given. Section |3] discusses the Monte Carlo 
(MC) simulation used. Section |4] describes the event selection and analysis methodology. The 
systematic errors are addressed in Section |5] and the results are presented in Section |6j 

2 The CMS Detector 

A detailed description of the CMS detector can be found in Ref. ||2T| . The subdetectors used 
for the present analysis are the inner tracker, consisting of silicon pixel and silicon strip layers, 
and the muon detectors. The inner tracker is immersed in a 3.8 T axial magnetic field. The 
pixel tracker consists of three barrel layers and two endcap disks at each barrel end. The strip 
tracker has 10 barrel layers and 12 end-cap disks at each barrel end. Muons are measured 
in gas-ionization detectors embedded in the steel return yokes. In the barrel, there is a drift 
tube system interspersed with resistive plate chambers (RPCs), and in the end-caps there is a 
cathode strip chamber system, also interspersed with RPCs. The first-level (LI) trigger used 
in this analysis is based on the muon system alone, while the high-level trigger (HLT) uses 
additional information from the inner tracker. 

The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
LHC beam collision point, the x axis pointing towards the center of the LHC ring and the z axis 
pointing along the counterclockwise beam direction. The polar angle 9 is measured from the 
positive z axis and the pseudorapidity is defined by = — In tan(0/2). The azimuthal angle (p 
is measured from the positive x axis in the plane perpendicular to the beam. 
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4 Data Selection and Analysis 



3 Monte Carlo Simulation 

The MC event generator PYTHIA 6.422 1221 is used (with MSEL=1) to compute efficiencies and 
kinematic distributions. PYTHIA and MC@NLO 3.4 11231 124| predictions are compared with the 
experimental results. The programs were run with their default parameter settings, except 
when mentioned otherwise. The PYTHIA event sample was simulated with the CTEQ6L1 [25^ 
PDF, a b-quark mass m\, = 4.8 GeV, and Peterson et al. fragmentation functions [26J for c and 
quarks with parameters = 0.05 and ey = 0.005. The underlying event is simulated with the 
D6T tune [27j. Pileup events were not included in the simulation and play a negligible role in 
the data sample used for this measurement. 

For comparison, additional event samples were generated where the evtgen [28J program was 
used to decay the b hadrons. Events generated by the PYTHIA program were passed through 
a detailed MC simulation of the CMS detector response based on GEANT4 [291. The MC@NLO 
package has a NLO matrix element calculation interfaced to the parton shower algorithms of 
the HERWIG |30| package. A b quark mass of nth = 4.75 GeV and the CTEQ6M PDF set |23 
were used. The events generated with MC@NLO are studied only at the generator level and are 
not passed through the detailed detector simulation. 

4 Data Selection and Analysis 

This analysis is based on data collected in 2010 when the collider and detector were fully oper- 
ational and fulfilled the following requirements: (1) Stable beam conditions, (2) stable magnetic 
field inside CMS at the nominal value, (3) operational LI and HLT, and (4) inner tracker and 
muon stations at their nominal high-voltage settings. The data sample used in this analysis 
corresponds to an integrated luminosity of £ = 85 ± 9nb^^ [31 1. 

The events of interest are selected by a very loose single-muon trigger path. The LI muon 
trigger makes no explicit requirement on the muon momentum transverse to the z axis, pj, 
although muons with pj < 3 GeV do not have sufficient momentum to be reconstructed in the 
barrel region of the muon system. 

In the HLT, a standalone muon reconstruction (with information from the muon detectors only) 
is seeded by the parameters of the LI muon candidate. If the standalone muon candidate has 
Pj > 3 GeV it serves as a seed in the global muon reconstruction, where a track in the inner 
tracker is linked to the standalone muon, and further selection requirements are applied on the 
transverse momentum {pj > 3 GeV) and the impact parameter with respect to the beam spot 
in the transverse plane (|do| < 2 cm). 

The offline event selection requires a reconstructed primary vertex with more than three tracks 
and at least one muon candidate with pj > 6 GeV and pseudorapidity < 2.1 that fulfills a 
tight muon selection similar to that in Ref. ||32| . The muon candidates are required to be recon- 
structed by two independent algorithms, one starting from segments in the muon chambers 
and one starting from inner-tracker information. The inner track must be measured with at 
least 10 hits in the inner tracker, two of which must be on pixel layers. The inner-track fit and 
the global muon fit (including all inner tracker and muon detector hits) are required to have a 
of less than 10 per degree of freedom and at least two muon segments matching the inner 
track must be found. Only muon candidates with transverse impact parameter with respect to 
the primary vertex \do\ < 2 mm and longitudinal impact parameter with respect to the primary 
vertex \dz\ < 1 cm are accepted. 

In events passing the trigger and event selections, all tracks including the muon are clustered 
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into track-jets l|33| by the anti-Zcx jet algorithm [34] with R = 0.5. The tracks are selected with the 
following requirements: 0.3 < px < 500 GeV, |zo | < 2 cm, and hits in at least 2 (5) layers of the 
pixel (pixel and strip) detector. Only jets containing a muon are accepted as b-jet candidates. 

The jet direction and jet energy E are calculated by summing the four-momenta of all tracks 
in the jet except the muon. The pion mass hypothesis is assumed for calculating the energy 
associated with a track. The jet is required to contain at least one track and to have a transverse 
energy Ej = E sin 0jet of at least 1 GeV, where 0jet is the polar angle of the jet direction. 

The efficiency for identifying b jets is determined in MC simulation for events in which the 
muon from a b-hadron decay falls into the kinematic region of this measurement. The efficiency 
for finding a jet containing the muon rises with the muon pj from 74% at 6 GeV to almost 
100% for events containing a muon with pj > 20 GeV. The fraction of events in which the 
reconstructed jet containing the muon is not matched to the b jet at the generator level is smaller 
than 7% in the lowest muon transverse momentum bin and asymptotically reaches a value of 
2% at large pj. 

From the momenta of the selected muon (p^,) and the associated track-jet (py), the relative 
transverse momentum of the muon with respect to its track-jet is calculated as p'^' = |p^ x 
Pj\/\Pi\- 



x10^ 




muon pf [GeV] 



Figure 1: Distribution of muon transverse momentum p"^ with respect to the closest track-jet 
in data and results of the maximum likelihood fit. The black full circles correspond to the data 
distribution, while the black line is the result of the fitting procedure. The red dashed and the 
blue dotted line are the simulated b and cudsg distributions, respectively. 



A total of 157783 data events pass the selection. If an event contains more than one muon of 
either charge, only the muon with the largest transverse momentum pj is kept. This affects 
0.5% of all data events. 
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5 Systematic Uncertainties 



4.1 Fitting Procedure 

A fit to the observed p'^f' spectrum, based on distributions obtained from simulation (signal 
and cc) and data (the remaining background), is used to determine the fraction of signal events 
among all events passing the event selection. A binned log-likelihood fit is performed, which 
takes into account the finite size of the MC simulated sample [35] . 

The distributions used in the fitting algorithm are determined separately for the full sample and 
for each bin in muon transverse momentum and pseudorapidity. Since the shape of the p'J^ dis- 
tribution in c and light-quark/ gluon (udsg) events cannot be distinguished by the fit, the two 
background components are combined and a fit discriminating the signal component against a 
single background component is implemented. The udsg background is dominated by hadrons 
misidentified as muons (mainly in-flight decays) and is determined in data. Hadrons satisfy- 
ing all muon track selection criteria (without muon detector requirements) are weighted by the 
misidentification probability and used instead of muons to determine The misidentifica- 
tion probability has been measured in data 1,36 J . The c background is determined from MC 
simulation. Muons from sources other than b, c and udsg events are neglected. The largest 
contribution to the muon event sample from these sources is expected in the highest pj bin 
(3%, from W decays). 

The result of the fit in the full sample is displayed in Fig. [Tj Extensive tests to validate the 
fitting procedure were performed [37] with repeated fits of MC pseudo-experiments obtained 
by appropriate random variations. A satisfactory performance of the fit was observed: the fit 
result does not show a significant bias and the errors are properly calculated by the fitter. The 
stability of the fit was proven by repeated fits with varied binning. The signal fractions have 
also been determined with particle flow jets |38| and with a fit to the muon impact parame- 
ter distribution. The results are cor\sistent with the fit using track-jets within the systematic 
uncertainty. 

5 Systematic Uncertainties 

The systematic uncertainties of this analysis are dominated by the shapes of the distribu- 
tions used in the fitting procedure. 

The signal p'^' distribution is validated with data through a control sample enriched in b de- 
cays. Selecting muons with a large impact parameter significance of |do I / ^dg > 12, where do is 
the uncertainty of the impact parameter measurement, results in an event sample with an ex- 
pected b fraction of about 85%. Small adjustments of the shape of the distributions by rescaling 
p^^^ improve the agreement between data and simulation in the b-enriched region and in the 
full sample. They result in variations of the measured cross section of up to 21% that are taken 
as a systematic uncertainty. 

The background consists of contributions from cc events and from light-quark and gluon events, 
where a hadron is misidentified as a muon. Both contributions are similar in shape and mag- 
nitude. The c fraction of the background is expected to rise with increasing muon pj. The 
fit does not separately determine the c and udsg content of the sample. Two effects can in- 
troduce a systematic uncertainty. (1) The udsg distribution determined from data could be 
biased. Using the PYTHlA-derived udsg background introduces a difference to the reference fit 
of 2-14%, depending on the muon transverse momentum and pseudorapidity bin. (2) If the c 
fraction of the non-b background in the data was different from the value used in combining 
the backgrounds, the fitted b fraction could change. The MC-simulation predicts a c fraction of 
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Table 1: Summary of systematic cross section uncertainties. The systematic uncertainty can 
vary depending on the muon transverse momentum and pseudorapidity as indicated by the 
range. 

source cross section uncertainty (%) 

Trigger efficiency 5 

Muon reconstruction efficiency 3 

Hadron tracking efficiency 2 

b p^f - shape uncertainty < 21 

Background p'J^ shape uncertainty 2-14 

Background composition 3-6 

Production mechanism 2-5 

Fragmentation 1-4 

Decay 3 

Underlying event 10 

Luminosity 11 



50-70% in the non-b background depending on the muon transverse momentum. This fraction 
depends on the modeling of charm semileptonic decays and on the hadron misidentificaton 
probability. Varying the c vs. udsg fraction by ±20% leads to a systematic uncertainty of 3-6%. 

The muon trigger efficiency has been determined from data with an uncertainty of 5% using 
independent triggers. The muon reconstruction efficiency is known to a precision of 3%. The 
tracking efficiency for hadrons is known with a precision of about 4% [39], which induces a 
systematic uncertainty of 2% on the number of events passing the event selection. 

In PYTHIA, the production of a bb pair can be separated into flavor creation (19% of the selected 
events), flavor excitation (56%), and gluon splitting (25%). The event selection efficiencies are 
71%, 72%, and 76%, respectively. Reweighting the events from the different production pro- 
cesses to reflect the difference between pythia and HERWIG leads to a systematic uncertainty 
of 2-5%, depending on the muon transverse momentum. The uncertainty of the b quark frag- 
mentation is studied by varying the parameter between 0.003 and 0.010, which results in 
a systematic uncertainty of 1-4% on the reconstruction efficiency. A sample generated with 
EVTGEN is used to investigate the uncertainty in modeling the b-hadron decay properties. A 
systematic uncertainty of 3% is found. Varying the fraction of prompt b — ?> decays with 
respect to b — > c — > decays within its uncertainty [40J changes the measured cross section 
by 1%. Neither the muon trigger efficiency nor the track-jet finding is affected significantly 
by the variation of the fragmentation and decay parameters. The track-jet reconstruction can 
be affected by the underlying event. Using simulated event samples with different MC tunes 
(D6T [27], Pro-Q20 ll4T]i , and CW [42]) for the efficiency and acceptance calculation changes the 
cross section of the order of 10 %. At the present stage of the CMS experiment, the integrated 
luminosity recorded is known with an accuracy of 11% [|31| . 

Table [l] summarizes the systematic uncertainties. 
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7 Conclusions 



6 Results 

The inclusive production cross section for b quarks decaying into muons is calculated as 

where N\, is the number of selected b events in data. No distinction is made betw^een posi- 
tive and negative muons; Nb includes the process pp ^b + X — > fi + X'. The efficiency e 
includes the trigger efficiency (88 ± 5)%, the muon reconstruction efficiency (94 ± 3)%, and 
the efficiency for associating a track-jet to the reconstructed muon, (77 ± 8)%. 

The result of the inclusive production cross section for b quarks decaying into muons within 
the kinematic range pj > 6GeV and \t]f'\ < 2.1 is 

cr = 1.32 ± 0.01 (stat) ± 0.30(syst) ± 0.15(lumi) }ih, 

where the first uncertainty is statistical, the second is systematic, and the third is associated 
with the estimation of the integrated luminosity. For comparison, the inclusive b-quark pro- 
duction cross section predicted by MC@NLO is 

(7-MC@NLO = 0.84+°:3^ (scale) ± 0.08(mb) ± 0.04(pdf) /Yb, 

where the first uncertainty is due to variations in the QCD scale, the second to the b-quark 
mass, and the third to the parton distribution function. The value of the scale uncertainty is 
obtained by varying the QCD renormalization and factorization scales as described in Ref. [Zl. 
The b-quark mass was varied between 4.5 GeV and 5.0 GeV and the uncertainty induced by the 
parton distribution function was evaluated using the eigenvector sets as described in Ref. II25I . 
The PYTHIA prediction using the parameters described in Section|3]is 1.8 }ib. 

The differential cross section is calculated from 

K 

bin,- ^^'A^'" 

where x stands for the muon transverse momentum or the muon pseudorapidity, and Ax' de- 
notes the width of bin i. The number of selected b events in data and the efficiency e, are 
determined separately for each bin. 

The results of the differential b-quark production cross section as a function of the muon trans- 
verse momentum and pseudorapidity are shown in Fig.|2]and summarized in Table |2] The data 
lie between the PYTHIA and the MC@NLO predictions. The observed shapes of the kinematic 
distributions are described reasonably well by both programs. The integral of the differential 
cross section is consistent with the cross section determined from the full sample. 

7 Conclusions 

A measurement of the inclusive b-hadron production cross section in the central rapidity region 
in proton-proton collisions at a/s = 7TeV has been performed. The measurement is based 
on a data sample corresponding to an integrated luminosity of 85 nb^^ recorded by the CMS 
experiment during the first months of data taking in 2010 with a low-threshold single-muon 
trigger. 
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Figure 2: Differential cross section (left) j^{pp ^h + X^]d + X', < 2.1), and (right) 

dp J 

1§F (PP —>h + X—>}i + X',pj>6 GeV) . The two possible muon charges are not distinguished 
and the process pp — > b + X — > + X' is included. The black points are the CMS measure- 
ments. Vertical error bars showing the statistical error are smaller than the point size in most 
bins, the horizontal bars indicate the bin width. The yellow band shows the quadratic sum 
of statistical and systematic uncertainties. The systematic uncertainty (11%) of the luminosity 
measurement is not included. The solid blue line shows the MC@NLO result and the dashed 
blue lines illustrate the theoretical uncertainty as described in the text. The solid red line with 
dots shows the PYTHIA result. 



The result for the total inclusive production cross section of b hadrons decaying into muons 
within the visible kinematic range is 

(7(pp ^ b + X ^ ^( + X') = 1.32 ± 0.01 (stat) ± 0.30 (syst) ± 0.15(lumi) jib, 

where pj > 6 GeV, \rjf\ < 2.1. The measured cross section is approximately 1.6 times higher 
than the MC@NLO prediction, but the difference is less than the theoretical and experimental 
uncertainties. Differential cross sections have been measured as a function of muon trans- 
verse momentum and pseudorapidity. The observed shapes are reasonably well described by 
MC@NLO. A similar pattern was recently found by this collaboration in the measurement of b 
production using fully reconstructed meson decays HI. 
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7 Conclusions 



Table 2: Differential cross sections da I dpj for \r]f'\ < 2.1 in bins of muon transverse momen- 
tum and da/ drj^' for pj>6 GeV in bins of muon pseudorapidity. The number of b events 
(Nb, including b events) determined by the fit, the efficiency (e) of the online and offline event 
selection, and the differential cross section, together with its relative statistical and systematic 
uncertainties, are given. A common uncertainty on the luminosity of 11% is not included. 
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